ANL-HEP-PR-11-67 
NUHEP-TH/ 11-24 



O 

(N 

> 
O 

o 

m 



o 



X 



A subtraction scheme for NNLO computations 

Radja BoughczaQ 

High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439, USA 

Kirill Mclniko\|l| 

Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 

Frank PctricllcH 

Department of Physics & Astronomy, Northwestern University, Evanston, IL 60208, USA 
High Energy Physics Division, Argonne National Laboratory, Argonne, IL 60439, USA 



We use the known soft and coUinear limits of tree- and one-loop scattering amplitudes - computed 
over a decade ago - to explicitly construct a subtraction scheme for next-to-next-to-leading order 
(NNLO) computations. Our approach combines partitioning of the final-state phase space together 
with the technique of sector decomposition, following recent suggestions in Ref. [l[. We apply this 
scheme to a toy example: the NNLO QED corrections to the decay of the Z boson to a pair of 
massless leptons. We argue that the main features of this subtraction scheme remain valid for 
computations of processes of arbitrary complexity with NNLO accuracy. 



I. INTRODUCTION 

Asymptotic freedom and factorization of short- and 
long-distance effects in the hard scattering of hadrons to- 
gether allow us to employ perturbative computations in 
QCD to extract information about the physics of parton- 
parton interactions. Perturbative QCD computations are 
organized by "loop order" in an expansion in the strong 
coupling constant. At both leading- and next-to- leading 
order in Us, the conceptual framevifork for perturbative 
QCD computations is well-established. While this does 
not necessarily make such calculations easy, or even fea- 
sible in some cases, the existence of a general framework 
is needed to claim a full understanding of the structure 
of perturbative QCD at these orders. 

The crucial element of such an understanding is the 
concept of infrared and collinear safety. It leads to com- 
putational algorithms Q that permit the calculation of 
arbitrary hadron-coUider observables to next-to-leading 
order (NLO) accuracy. These methods rely upon approx- 
imating matrix elements containing an additional mass- 
less particle by appropriate limits in singular regions of 
phase space. These limits are simple enough to be in- 
tegrated over the unresolved phase space of the radiated 
particle, without any reference to a particular observable. 
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Unfortunately, the situation is much more confusing at 
next-to-next-to-leading order (NNLO) and beyond. To 
illustrate this point, we note that no NNLO computa- 
tion of a cross section for a 2 2 scattering process 
with strongly interacting particles exists, in spite of the 
fact that a large number of two-loop virtual amplitudes 
[3-01 and all singular limits of tree- and one-loop ampli- 
tudes have been known for over ten years. The 
reason is that a working algorithm that combines these 
ingredients to obtain physical cross sections has not been 
formulated. Consequently, a significant number of ex- 
isting fully differential NNLO computations jlln25j were 
performed using unorthodox approaches, that are only 
remotely related to mainstream NLO subtractions ideas 
considered generalizable to higher orders . 

An alternative NLO subtraction scheme due to Frix- 
ione, Kunszt and Signer (FKS) [1^ , was not the scheme 
of choice when NLO computations were initially be- 
ing undertaken, but has received renewed interest re- 
cently [27|. The main idea of the FKS subtraction is to 
partition the phase-space such that in every sector only 
one definite external particle i can become soft, or two 
definite external particles i and j can become collinear to 
each other. If such a partitioning exists, it is clear that 
all singularities in a given sector are easily extracted if 
the unresolved phase space is parameterized in terms of 
the energy of particle i and the relative angle between 
directions of particles i and j. 

In the case of NNLO computations the "elementary 
building block" is the double-unresolved phase space, 
where two given particles can become soft or collinear 
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to a third particle. Extraction of singularities in the 
triple collinear limit is non-trivial, but can be accom- 
plished by applying the technique of sector decomposi- 
tion [ll|, I2I, [231 to the three-parton unresolved phase- 
space. However, early applications of sector decomposi- 
tion at NNLO [ill . [isl Il6l | did not perform an initial par- 
tition of the phase space to separate collinear singulari- 
ties, and instead attempted to find a suitable phase-space 
parameterization for an entire process at NNLO. The 
complexity of this endeavor slowed down the progress 
of NNLO computations for hadron-collider observables 
using sector decomposition after initial success with such 
computations for the hadro-production of the Higgs and 
electroweak bosons. 

It was recently pointed out by Czakon [l], [3^ that by 
combining the idea of phase-space partitioning from FKS 
with the idea that sector decomposition can be applied to 
real emission integrals [ill [29| , a powerful framework for 
NNLO computations is obtained. The purpose of this pa- 
per is to elaborate on this observation and show explicitly 
how this framework can be used to obtain physics results. 
In addition, we show how to handle collinear singularities 
in the final state, and discuss in detail the computation 
of real-virtual corrections and the scheme dependence of 
the results. For the sake of simplicity, we study two-loop 
QED corrections to the decay rate of the Z-boson to an 
electron-positron pair. Dealing with QED corrections of- 
fers significant simplifications, yet is far from trivial, and 
is a good place in which to develop and test ideas. 

In the following sections we explain in detail our ap- 
proach. In Section II, we discuss the extraction of sin- 
gularities from double-real radiation corrections. In Sec- 
tion III, we explain how one-loop corrections to the one- 
photon real-emission process are treated. In Section IV 
we describe the dependence of the results on the choice 
of the regularization scheme. In Section V we present 
our conclusions. 



We begin by describing the collinear partitioning. We 
introduce four functions 



A± = 1 



(1) 



where rii is a unit vector in the direction of the pho- 
ton i and are the unit vectors in the directions of 
the positively- and negatively-charged leptons. For each 
photon, we introduce a partition of unity 



1 = 



a; 



a; 



A+ + A^ A+ + A: 



= Pi + Pi : 



(2) 



and obtain 



1 - m 



+ Pi)^ Ptpt + P1P2 + P1P2 + Pi P2- (3) 



We introduce pfp'j = S- 
equation as 



1 



and re-write the previous 



(4) 



Each of the contributions on the right-hand side of 
Eq. (|4|) defines a primary sector whose phase space we 
must parameterize separately. The superscripts of each 
5 indicate the potentially singular collinear directions for 
each of the two photons. For example, in the primary sec- 
tor labeled by S^2~ both photons can become collinear to 
the electron. In the sector labeled by 6^2~ o^ily pho- 
ton can be collinear to the electron, while the other can 
become collinear to the positron. Using Eq. Q, we de- 
compose the phase-space as 



dLipSg_,_ 



e— 7172 



dLips°+ 



e — 7172 ' 



(5) 



a,b=± 



where 



II. DOUBLE REAL RADIATION 

We consider the decay of a Z-boson to an 6+6^77 fi- 
nal state, Z{pz) e+(p+)-|-e"(p_) + 7(pi)-f 7(^2)- We 
must design a strategy to integrate the squared matrix 
element for this process over the phase space of the final- 
state particles. This is non-trivial because such an inte- 
gration eventually leads to phase-space regions that con- 
tain collinear^ and soft singularities. The main problem 
is not that such singularities exist, but that in different 
regions of phase space, different subsets of particles lead 
to singularities in the matrix elements. As a first step, 
we partition the phase-space such that we know which 
final-state particles can develop singularities in each par- 
tition. 



dLips"'' 



e+e-7i7: 
ds-d 



- /[dp_][dp+][dpi][dp2] 



(6) 



X {2nf5\pz -P--P+-PI- P2)6f, 



^ We treat all final state particles as massless. 



For convenience, we introduce the short-hand notation 
[dp] = d'^-V/(2po(27r)''-i), where d = 4 - 2e is the di- 
mensionality of space-time. The overall 1/2! is the sym- 
metry factor for the two final-state photons. 

Out of the four primary sectors, two sectors contain 
triple-collinear singularities where both photon momenta 
are collinear to either the electron or positron momen- 
tum. The other two sectors contain doublc-coUinear 
singularities, where the momentum of one photon is 
collinear to the electron momentum while the other pho- 
ton is collinear to the positron. It is sufficient to under- 
stand one triple-collinear and one double-collinear par- 
tition. The remaining primary sectors are obtained by 
a simple re-labeling of the final-state particles. When 
discussing phase-space parametcrizations in the relevant 
sectors, we make use of the fact that the electron and 
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positron cannot develop soft singularities, and that in 
the process Z e"'"e~7i72, the kinematic configuration 
where the electron and positron momenta are coUinear 
is non-singular. If such singularities could occur, only 
an additional partitioning would be required to handle 
them. 



A. The triple collinear sector 

1. Phase-space parameterization 

We consider the S^2~ primary sector, where the pho- 
tons and the electron can develop collinear singularities. 
We must also consider soft singularities that appear when 
the energies of one or both photons vanish. Our discus- 
sion of the phase-space parameterization closely follows 
Refs. [ll,|33l. 

We first explain how the energies of the two photons 
are parameterized. We denote the sum of the four- 
momenta of the electron and positron by Q = p+ + p- . 
Momentum conservation implies pz ~ Pi ^ P2 ^ Q and 
< < to|. We write to| - ^ Am|, < A < 1. 
Squaring the momentum conservation equation, we find 

m| - 2mz{Ei + E2) + 2EiEi{l - m ■ ns) = (3^ (7) 

where E12 are the energies of the two photons and rii 2 
are three-dimensional unit vectors along their momenta. 
We parameterize the photon energies by Ei = S^iinz/I, 
and the relative angle between them by r/12 — {\ — ni ■ 
n2)/2. Solving for ^1 or ^2 in Eq. ([7]) then yields 



6 = 



A-e2 

1 - 6'n2 



or ^2 



A -6 

1 - 6'7l2 ' 



(8) 



We can remove the symmetry factor in Eq. ([6|) by re- 
quiring that El > E2 and by using the fact that the 
matrix element is symmetric under the interchange of 71 
and 72. We obtain 



dLips, 



'e+e— 7i72 



{27rf / [dp.][dp+][dpi][dp2K, 



where 

1-a 



^max(Ci) = min 



1, 



6(1 - 6'7i2) 



(9) 



(10) 



We decompose the four-particle phase-space into "regu- 
lar" and "singular" phase-spaces: 



dLipSg+e_^^^2 = dLips^eg x dLips^; 



(11) 



where 



dLipS,::„- = [dpi] [dp2] 5^2~ ^(a^max(a) " 6), (12) 



and 



dLips,,^ = [dp_][dp+](27r)'^<5(^)(Q (13) 



We begin with a discussion of the singular phase-space. 
We note that since the Z-boson decays at rest and there 
arc only four particles in the final state, we can choose 
the momenta of any three particles to be in the four- 
dimensional space, without any (d — 4)-dimensional com- 
ponents. The three- momentum of the fourth particle is 
determined by momentum conservation, and is also in the 
four-dimensional space. To have a simple parametriza- 
tion, we choose the direction of the electron momentum 
to be the z-axis. Then, p_ = i?_(l,n_), n_ = (0,0,1) 
and pi,2 = £^1,2 (1, JT-i), where rii = (sin6'i, 0, cos^i) and 
722 = (sin 6*2 cos sin 6*2 sin cos 02 )■ We also introduce 
the following notation for the scalar product of n_ with 
"1,2: 



^71,2 



1 - n_ • ni,2 



(14) 



We then find the scalar products 

2pi • P- = 2E_mz£,ini, 2pi ■ p2 = mz^i£,2Vi2- (15) 

We now write the parametrization of the singular phase- 
space using the angular variables just introduced: 



,T. — ^^12 ^(Ci -6) fmz\ 



2d-4 



X dn^f-'^dn^^^'^ d6 cj-'^ d6 e2~^' 



'2 (16) 
X d77i [771(1 - 77i)]"M772 [772(1 - ?72)]~' 

X 



d cos f [l — cos^ ip 



2 „\-l/2-£ 



The goal is to rewrite the integration over cos ip in a way 
that makes the factorization of singularities manifest. To 
this end, we introduce the variable k: 



(l-cos(6>i -6'2))(l + cos(^) 
2 (1 — cos(0i — 6'2) + (1 — cos Lp) sin 6*1 sin 02 ) 



(17) 



Because cos (p can be used to parameterize 7^12 , we will 
need a relationship between 7712 and k. It can be easily 
derived by solving the above equation for cos p and then 
using the solution in the expression for 7712. We find 



?7l2 



where 

N{rii,r]2, k) V2 ~ 2771772 



2(1 - 2k)v/77i'72(1 -??i)(l -r/2). 



(18) 



(19) 



Finally, we need the Jacobian for the ip ^ n variable 
transformation, and a simple expression for sin^ ip = 1 — 
cos^ ip. The relevant equations are 



d cos tp 277^2 



dn 



1 — cos ip 



2 Ak{1 - K)ril2 



(20) 
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We can now change variables </? — s- k in Eq. (IT6l) for the 
singular phase-space. We find 



dLipSsi 



2d-4 



2d-2 



smg 24+2^ (27r) 

X d?7i [?7i(l - ?/i)]"'d?72 [772(1 - 772)]"' 



(21) 



X dK (k(1 — k)) 



l-2e 
-1/2-e ^12 



I??! - ml 



As we will see later, Eq. (PT|) gives us the singular phase 
space in a form that is convenient for the extraction of 
singularities. 

Next, we discuss the regular phase space. We write it 



dLips_ = [dp.][dp+]{27:fS^''HQ-p+-p^) 



dfi 



(d-i) 



E 



l-2e 



(22) 



2(27r)''-2 2(0o - Q • n_ 



where = Q'^ /[2{Qo — Q-n^)] is the electron energy. It 
is important to understand which elements of the calcu- 
lation can be simplified by setting the number of space- 
time dimensions to four, d — !■ 4. In the context of the 
phase-space discussion, we factor out the leading order 
phase-space for Z — >■ e+e" in Eq. ([22]) and treat is as 
four-dimensional. Everything else in Eq. ([22]) is treated 
with exact e-dependence. The leading-order phase-space 
that we use in what follows reads 



dLips 



df]('^-i) /mz\-2e dcos 
~ ^ 327r2 



8(27r)'^-2 V 2 
The regular phase-space becomes 

2_E_ 

dLips = dLips^. 



(23) 



{Qo-Q- n- 



2E^ 
mz 



(24) 

where dLips2_j.g+g- is taken in four dimensions, as in 
Eq. 

With the explicit parametrization of the phase space 
at hand, we are ready to discuss how the momenta of the 
final-state particles are generated. We follow the simple 
procedure described below: 

1. first, we use Ci,2, ?7i,2, to generate the energies and 
momenta of the two photons; 

2. second, we calculate Q = pz ^ Pi ^ P2] 

3. third, we obtain = Q^/{2{Qq — Q ■ «-)), as- 
suming n_ is along the z-axis; 

4. finally, the four-momentum of an electron is taken 
as p_ = -E- (1, 0, 0, 1), and the four-momentum of 
the positron is calculated by = pz —pi —p2 —P- ■ 



The remaining issue is the extraction of singularities 
from the matrix element of Z — ^ e"'"e~7i72 in the S^2~ 
sector. The matrix element is singular if either photon 
is soft, and also if either photon momentum is collinear 
to the electron momentum. By analyzing the potentially 
singular denominators, it is straightforward to find that 
there are three sectors to consider. These sectors are 
identified by changes of variables (5, 77, k) — >• {xi^i,,^} 
that we make in order to factor out all singularities from 
the matrix element. The three sectors are^: 



1. , where ^1 = xi, ^2 

772 = XiXs, K = X5] 



XmaKX2Xl, ?7l = X3, 



2. S2 , where = xi, ^2 = Xmi,^X2X4Xi, rji = xsXa, 

r]2 =X3, K = X5; 

3. 53"", where = xi, ^2 x,na.^X2Xi, 771 = 2:2X3X4, 

772 = X3, K = X5. 

For each of the sectors Si we must express the phase-space 
through the new variables and find the singular limits of 
the amplitudes. We illustrate how this is accomplished 
for the sector . The remaining two sectors are han- 
dled in a similar fashion. 



2. The sector 

For the sector S^^^ we write the phase-space in the 
following form 

dLipSg^^ = dLips"^^ ^xfx2x'^X4m^S^2'] 1 



dLips V = dNorm PS„ (PS)" 

O 1 

dxi dx2 dx3 dx4 



dn 



(25) 



xl+^' xl+^' xl+^' xl+' 7r(«;(l - A^))i/2 ' 

where 

dNorm = dLips^_,,,- liilil!^^^ , 

Bf^ = l-^e'-2C3e' + ^e^ + 0{e% 
and the normalization factors are given by 

2;max(l - X4) 2E^ 



(26) 



PS„ 



N{x3,X4, X5) {Qo Q n^) 
n 2 



PS = 16 



(1 - X4) 



N{x3,X4,X;i) 



{I - X3){1 - X4X3) (27) 



X [k(1 - K)]2;^ax ( 



Note that the normalization of the various pieces is cho- 
sen so that the TTi-^-dependent factor in dNorm can be 



Wc use the notation Xmax for the function 5max(a:i) 
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factored out entirely. This will be true for both the two- 
loop virtual and real-virtual contributions. The factor 
of m| in the square brackets in Eq. (f25|) is present to 
make the product of the square bracket and the ampli- 
tude squared for Z — > e+e~77 dimensionless. Note also 
that dNorm in Eq. ([26| is the same for all sectors, and 
that we will be using it in equations for other sectors be- 
low. The function iV(x'3, X4, 2:5) appears in the equation 
for 7712 after it is expressed through the new variables. It 
reads 



2(1 - 2x^)^x^{\ - xz){\ - xzx^. 



(28) 



The momenta of all the particles arc written through x- 
variables in the following way: 



?/l2 = 



a;3(l - X4y 



N{x3,Xi,X5)' 

cos 6*1 = 1 — 2a;3, cos 6*2 = 1 — 2x3X4. 



(29) 



Angles 61^2 are the polar angles for the two photons. We 
take sin 6*1^2 to be positive-definite. We choose photon 71 
to be in the x — z plane, so that (^1 = 0. The azimuthal 
angle of the photon ip2 is calculated to be 



sm (p2 



\/4:X5{l - X5){1 - Xj) 
N{X3,X4,X5) 



, COS tp2 ~ Xy 1 — sin tp2, 
(30) 

where A = sgn((l — cos0icos02) — 27712). Finally, the 
momenta of the two photons are given by 



Pi = Ei{l, sin^i, 0, cos^i), 

P2 ~ i?2(l, sin 02 cos (/?, sin 02 siniy9, cos Q2) 



(31) 



The momentum of the electron is p_ = ii^_(l, 0, 0, 1), 
and the positron four-momentum is obtained from mo- 
mentum conservation. 

With the parameterization of the phase-space at hand, 
we can discuss extraction of singularities from the matrix 
element. To this end, the factor in square brackets in 
Eq. (pS]) is combined with the amplitude squared. This 
should give a finite expression in all singular limits. We 
introduce the regular function 

FliXXi^x..^}) = \x\x\x\x4m\b'{2\ IMz^e+e-fll'^ , 

(32) 

A calculation of Z — > e^e 77 contribution to Z-dccay 
rate involves integration of the function Fi over the 
phase-space in Eq. ([25]) : 



y d LipS ^^Fi{xi,X2,X3,X4,X5). 



(33) 



Because of the structure of dLips in Eq. the in- 

01 

tegrand is singular if one of the integration variables Xi , 



i = 1, ..4 vanishes. Such singularities can be extracted 
by writing 



1 



6ix)- 



' 1 ' 




Ina;^ 




— ne 




_Xi_ 


+ 


Xi 



-0(e2) (34) 



for all of the singular variables. The plus-distributions 
are defined in a standard way: 



1 i 
J dx — f{x) ~ j dx 



/(•^) - /(O) 



(35) 



Resolving the plus-distributions as in Eq. psp . we find 
that the integration in Eq. psp requires knowing Fi when 
one or more of its arguments vanishes. The calculation 
of these limits requires care. Although the function Fi 
is regular everywhere, the matrix element Aiz^e+e-^y^ 
is singular when certain a;- variables vanish. The physical 
limits that correspond to a set of particular x-variablcs 
vanishing can be deduced from the expressions for the 
momenta given in Eqs. p9l[3T|) . We describe the relevant 
limits below. 

• If xi = 0, the energies of both photons vanish. This 
corresponds to the double-soft limit. In the double- 
soft limit, the QED eikonal currents completely fac- 
torize and we obtain 



\M 



Z^e+e 77 1 



e^JiJ2\Mz- 



(36) 



where the square of the eikonal current for the pho- 
ton i reads 



2p- • V+ 



{P- ■Pi)iP+ - Pi) 



(37) 



The scalar products arc computed using the ex- 
plicit paramctrization of the momenta given in 
Eqs.(l2l|30l[3T]). We obtain for the function Fi 



^lUi=0 



16e^ 



(38) 



• If xi 7^ but X2 = 0, the photon 71 is hard and the 
photon 72 is soft. The matrix element becomes 



\M 



Z— >-e+e 7172 



^e\j2\Mz- 



•■e~^e 71 I 



(39) 



in this limit. Calculating the function Fi for xi 7^ 
0,X2 = 0, we obtain 



-^"112:2=0 



>4(p+ ■p-)xlx3 (1 - n+ • ni) 



X \Mz^e+e-^^f, 
2 

where A12 — Y[ (1 " ''^-i • ' 

i=l 



A12 



(40) 



6 



• Equation (|40p develops singularities when — >■ 0, 
in which case photon 71 becomes collinear to the 
electron. We do not show hclicity labels in what 
follows because helicity is conserved along massless 
fermion lines. In the collinear limit we find 



\M', 



2e^ 

Sic 



Pe^{e,z)\M 



(41) 



where the momentum of is given by the sum 
P- + pii Pe-yi^, z) is the e — > e + 7 splitting func- 
tion given in the Appendix, and sie = 2p_ ■ pi = 
2E^mzXiX3. Upon evaluating Fi in that limit, we 
obtain 



^lU2=0, 1:3=0 — 



We'^xi 



mzE-xl^^^Ai2 



PeT,(e,z) 



(42) 



X \M2 



The fraction of energy carried away by the electron 
in the e — > e + 7 splitting is expressed through the 
variable z = 1/(1 + mzXi/{2E_)). 

• We next consider the X4 = limit, which corre- 
sponds to the photon momentum p2 being collinear 
to the electron momentum p_ . We calculate 



FI 



X4—0 



e'^XiX2X3Si2 fnz 



E^x^ 



(43) 



X \M 



iv— >e+ e 71 I 



where the momentum of e 
1/ (1 -I- mzx^^^X2Xi/{2E^)). 



is p_ -f p2 and z 



• Finally, we consider X3 = 0. This limit corresponds 
to the triple collinear limit, when the momenta of 
photons 71 and 72 are parallel to the electron mo- 
mentum p_ . In the triple collinear limit the matrix 
element factorizes as 



X P. 



2e^ 

Sl2e 



(44) 



,{e,Zf,,Zi,Z2)\Mz^e+e- 



where the momentum of is p- -\- pi + P2 and 
the triple splitting function ^67172 can be found in 
the Appendix. The energy fractions in this case are 
given by 



1 



mz 
2E- 



z\ 



mzxi 
2P_ 



XI (1 



Z2 



^max'^2) 1 
mzXiX2Xma.x 

2E^ 



(45) 



To compute Pi at X3 = 0, we introduce the notation 
si2e ~ 2E^mzXiX3di2e + 0{xl), 

dl2e = 1 + X2Xma.^Xi 

mz XiX2Xnie,x{l - ^4)^ 



(46) 



2P_ 



N{0,X4,X5) 



Using this notation, the function Pi is easy to write 
down. As an example, we present an explicit result 
for Pi(xi, a;2, 0, X4, X5) in the CDR regularization 
scheme: 

4 r 2 



2t 



2 2 

+i^2(e, Ze, Z2, H F2{e, Ze, Zi,Z2) 

2 2 2 

X-]X2 -r^ / \ ■.nX-\XnX4. 

+ -r-^ P2{e,Ze,Z2,zi)-{l-e)^ ^ ^ 



di2eXn 



X Xn'La-xX2X4 + 



dl2e 

+ 2e{l-ef^ 



(47) 



The functions Pi and P2 are presented in the Ap- 
pendix. 

The above formulae describe all the QED singular lim- 
its in the sector Si~. They can be implemented in a 
computer code in a straightforward way. 



B. The double collinear sector 



We next consider the i5{^^ primary sector. The singu- 
larities in this sector arise when the photon 71 is collinear 
to the electron, and the photon 72 is collinear to the 
positron. Soft singularities for both photons are possi- 
ble. 



1. The phase-space parametrization 

To develop a suitable description of the four-particle 
phase space in this partition, we make use of the momen- 
tum parametrization developed by Catani and Seymour 
for next-to-leading order calculations. Specifically, we 
adopt the momentum mapping for final-final dipoles, us- 
ing the language of Ref. l2| . We will also use the nomen- 
clature employed in Ref. [2| to describe different particles 
contributing to the dipoles, to make clear the connection 
to the discussion in that reference. 

We parameterize the phase-space for Z — >■ 6^6^7172 
in two steps. In the first step, we treat photon 71 as 
"emitted" , the electron as "emitter" and the positron as 
"spectator" . The 3^2 momentum mapping in this case 
is given in Ref. It is determined by the momentum 
conservation equation 

-1i+P-+p+=Pi-+p+, (48) 

and the relations between old and new momenta: 

pi = zipi^ +yi{l- zi)p+ + pi,x, 

P- = (1 - ^i)Pi- +2/i2:iP+ (49) 

P+ = (1 - Vi)p+- 
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The momenta pi_ and p+ are light-like: _ = = 0. 
The momentum pi ± is orthogonal to both of them. For 
the momentum parametrization in Eq. (|49| . the phase- 
space reads 



dLips(p_,p+,Pi,P2) = 



dLips(pi_,p+,p2) 



2! 



dyidzidil|j^i2 



4(27r) 



;2pi_ •_p+) 



(50) 



X {1 - y,)'-'^ yr (^i(l-^i))" 



In the second step, we apply a similar mapping for the 
momenta of the "reduced" reaction Z — ^ pi_ + p+ +P2, 
by considering 72 as "emitted", e"*" as "emitter" and ei_ 
as " spectator" . The momentum conservation equation 
becomes 

Pl-+P++P2=P+2+Pl-, (51) 

and the new momentum parameterizations read 

P2 = Z2P+2 +2/2(1 - Z2)Pi_ +P2,±, 

P+ = (1 - 22)^+2 + V2Z2P1- - P2.±, (52) 

Pi- = (1 - y2)pi-- 

Continuing with the parametrization of the phase-space 
shown in Eq. (|50p . we obtain 



dLips(p_,p+,pi,p2) = 



dLips(pi_,p2+) 



2! 



2^2-2e d^/ld2:ld^}|;lj2 d;/2dz2dl7y_^3 
X (^z) 4(2^)^-1 4(27r)'^-i 

X (l-2/2)'-'^(l- 22)^-^^(1 -2i)-^ 

We can express the momenta of all particles through 
2/1.27 •21.2- We choose Pi_,P2+ to be along the positive 
and negative z-axis respectively. The corresponding mo- 
menta read 

li„ = ^(1,0,0,1), (1,0, 0,-1). (54) 

We choose P2^± along the cc-axis and use Eq. ([5^ to ob- 
tain p+ and p2 ■ We find 

mz ( I — -. r — 

P2 = — [^2 +^2(1 - Z2),2^yZ2[l - 2:2)2/2, 

0,2/2(1 - Z2) - Z2), 

mz/ , (55) 

P+ = — (^(1 - ^2 +2/222), -2^22(1 - 22)2/2, 
0,2/222 - (1 - 22)). 

We can write momenta of the photon 71 and the electron 
using the above equations. The only semi-intricate step 
is the derivation of pi^±. We find 



l(2) 



(53) 



Pi,± = -7712^2122(1 - 2i)yi(l - 2/2)2/2"-+ 
+ mzV ^li^ - 21)2/1(1 - 2/2)(l - 22)n^, 



(56) 



where 
n 



— (cos (p, 0, 0, cos (f) , = (0, cos (p, sin 0) , 



Using Eq. (j56p . we derive the following expression for the 
energy of the photon 71 : 



2 



21(1 - 2/2) + 2/1(1 - 2l)(l - 22 + 2/222) 



(57) 



2a/?/i2i(1 - 21)222/2(1 - 2/2) COS(^ 



The energy of the photon 72 can be read off from Eq. (|55p . 

We now rewrite Eq. ((53)) by factoring out the Born 
phase-space and several other factors, similar to what 
has been done for the triple-collinear sector. From the 
momentum parametrization we see that, in addition to 
2/1,2,21,2, we need the azimuthal angle if to describe the 
phase space. There are no singularities associated with 
this angle in the double-coUinear sector, so we take it 
to be = 2Trx5, < < 1. In Eq. ([5^ we identify 
dLips(pi_,p2-i-) with the leading order phase-space and 
obtain 

dLips(p_,p+,pi,p2) = dNorm Psj„ PS'^^ 

_____ ^ I'^^j 

X d2/id2id2/2d22dx5y;^ "z^ '2/2 "22 '"^zi 
where dNorm is given in Eq. . The other factors read 

PS^ = i(l-2/i)(l- 2/2)^(1-22), 

PSgc„= 4(1 -2/1)^(1 -2/2)^^(1 -22)^ (59) 

X (1 — 2i)(l — COS^ if). 

There are two scalar products that can become singular 
in the coUinear limits: 



2p+ -^2 = (1 - 2/1)2/2™!, 

2p_ -pi = yi{l - j/2)(l - 22)to|. 



(60) 



Analyzing these scalar products and the expressions for 
the photon energies in Eqs. (j55l57p . we conclude that for 
the photon i = 1,2, the soft singularity corresponds to 
2/i ~ 0, Zi = 0, while the collinear singularity corresponds 
to yi = and 7^ 0. We also note that the apparent van- 
ishing of scalar products and photon energies atyi,Zi = 1 
also implies vanishing of the electron or positron energy, 
and for this reason does not lead to non-integrable sin- 
gularities. 

The two singular limits for each photon are factorized 
with the help of the two sectors 2/i < Zi and Zi < yi. Since 
there are two photons, we get four sectors altogether. We 
show below the changes of variables needed to completely 
factorize singularities in each sector: 

where 2/1 = xi, 21 = a;ia;2, 2/2 = 3:^3, 22 = 2:3X4, 

52"*", where 2/1 = a;2a;i, 21 = xi, 2/2 = 2:3, 22 = X3X4, 

where 2/1 = zi, 21 = X1X2, 2/2 = a;3a;4 22 = X3, 

54 ^, where 2/1 = X2X1, zi = .ti, 7/2 = a;3a;4, 22 = X3. 

(61) 
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Wc note that in sector Si~^, there are only soft singu- 
larities. In sectors 82^, , there is a soft singularity 
for one of the photons, and both soft and coUincar sin- 
gularities for the other. In sector S^^ there are soft 
and collinear singularities for both photons. The extrac- 
tion of all the limits in all sectors is similar to the triple 
collinear limit that we already discussed. For illustrative 
purposes, we only discuss the most difficult sector S^~^. 



2. Sector 



To discuss singular limits in sector S^^ it is convenient 
to introduce the following short-hand notation for the 
photon energies: 



El = —xiili, E2 = — a;3i22. 



(62) 



where 



fil = (1 - V2) + X2{1 - Zi){l ~ Z2+ 2/2^2) 



- 2^/2:2(1 - 2:1)2:22/2(1 - 2/2) C0S(/3, 
0,2 = I + X4^{1 - Z2) ■ 



(63) 



In this sector, singularities occur if any of the variables 
Xi, a;2, X3, 0:4 vanishes. To enable extraction of singulari- 
ties, we write the phase space as 



dLips, 



e^e 7172 



dNorm Ps,„ PS" 



dxi dx2 dx3 dx. 



l+2e 1+e l-|-2e 1 + 
1 2 3 4 

X [xiX2X3X4TO2'^f2^] 



(64) 



where dNorm is the same as in the triple collinear limit. 
The other factors read 



• If Xi = and X4 = 0, the photon 71 is soft and the 
photon 72 is collinear to the positron. The function 
F4 reads 



^4Ui=0,a;4=0 



32x36" 



(1 - 22)"l|f^lAi2 

X Pe7(e,^)l>^z^^+e-|^ 



(68) 



where 2 1/(1 + _E2/£'+) and the momentum 
is p+ +P2• 
• If Xi = 0, the photon 71 becomes soft and the func- 
tion i^4 reads 



8x3X4e^(p_ • p+)(l - n_ • n2) 
(1 - 2;2)(1 - Z2)mzE^VLxl\xi 



(69) 



X \M 



Z— ^e+e^72 I 



A new type of singular limit corresponds to X2 = 
and X4 — 0. This double-coUinear limit corresponds 
to photon 71 collinear to the electron and photon 72 
collinear to the positron. The function F4 evaluates 
to 



_ 16e"^xiX3 

-^4 2:2=0, 3:4=0 T7i 

?7l|(l - 22)Ai2 

X Pe7(e, 2l)^'e7(e, 22)|Xz->e-e- 



(70) 



where 21 = l/(l-f-£'i/-B_) and 22 = 1/(1+^^2/-^+). 
The momentum of e~ is the sum of the and 71 
momenta, and the momentum of is the sum of 
the e"*" and 72 momenta. 

The calculation of other limits proceeds along similar 
lines. The resulting expressions are again straightforward 
to implement in a computer code. 



PS» = i(l-2/i)(l- 2/2)^(1 -22), 

PS = 4(1-2/1)2(1 -2/2)^(1 -21) 
X (1 - 22)^(1 - COS^ 1^). 



(65) 



The finite function in this case is given by the product of 
the term in brackets in Eq. (j64p and the squared matrix 
element for Z — >■ 6+6^7172; 

^4({a;j=l,..5}) = XiX2X3X45f/|A/z^e+e-7i72p- (66) 

We now describe several of the singular limits in this 
sector. 

• The double-soft limit corresponds to Xi = and 
X3 = 0. In this limit, the function F4 evaluates to 



F4U,=o,.3=o = 20 o l-^^^e+e- r (67) 
m^"i"2 



III. VIRTUAL CORRECTIONS TO SINGLE 
PHOTON EMISSION 

In this Section we discuss the one-loop corrections to 
the single photon emission process Z — > 6+6^7. We 
first explain the partitioning of the phase-space and its 
parametrization. We then discuss how to compute sin- 
gular limits of scattering amplitudes. 

Consider the process Z — >■ e~e+7i. Singularities can 
arise if the photon is either soft or collinear to the electron 
or positron. To account for the collinear divergences, we 
partition the phase space as 



1 = + ^t. 



(71) 



where ~ pf , and are defined in Section |TT1 We 
write 



dLips 



e+e — i-i 



a=± 



dLips" 



(72) 
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where 



dLipSe+e-71 = J [dp-][dp+][dpi] 



(73) 



In what fohows, we consider the primary sector 6i , 
where a cohinear singularity can only occur when mo- 
menta of the electron and the photon become parallel. 
The other sector gives a symmetric contribution that 
can be analyzed identically. We use a parametrization 
that is similar to the two-photon case. For dLips~_|_g_^^ , 
we parameterize the photon energy as Ei — mzCi/'2 and 
the relative angle between the photon and the electron as 
cos 01 = 1 — 2?7i. The reference frame is fixed by requiring 
that the electron momentum is along the z-axis and that 
the photon momentum is in the x ~ z plane. Explicitly, 
we write 



p^^E^ (1,0,0,1), 
Pi = 



— (1, sin 01, 0, cos 6*1) . 



(74) 



The momentum of the positron is determined from mo- 
mentum conservation, = pz — P- — Pi- The en- 
ergy of the electron is found by first computing the 
momentum Q = pz — Pi and then calculating E- ~ 
Q'^ /[2{Qo — Q ■ rie)]. Finally, with this parametrization 
of the momenta, and borrowing notation that we already 
used when discussing the double-real emission, we write 
the phase space as 



dLips^_g+^ = Lips2_ 



(d-2) 



(2^) 



d-l 



ml'^'E- 



f2E^ 



-2e 



2(Qo Q n ) \ mz 
dCi di]i ,- 2 



(75) 



l + 2£ l + £ 
'/I 



We use Eq. ([75|) to construct a finite, integrable func- 
tion when it is combined with the squared matrix el- 
ement. To this end, we note that ^1 and r/i are al- 
ready suitable variables for the extraction of singularities, 
with ^1 controlling the soft limit and 771 controlling the 
coUinear limit. The squared matrix element in the real- 
virtual case is given by the interference of the tree and 
one-loop Z — >■ e~^e~j amplitudes. We employ Passarino- 
Veltman reduction to express the one-loop scattering am- 
plitude Z — !• e+e~7 in terms of one-loop integrals, and 
use the QCDloops program [3]| to compute master inte- 
grals. 

At first, it appears that we must simply repeat what 
we have done for the double-real emission corrections, 
extracting singularities of the matrix elements when ^1 
or 7/1 goes to zero. However, there is a subtlety here. 
One-loop amplitudes are not rational functions of ^1 and 
771, in contrast to their tree-level counterparts. To ame- 
liorate this problem, we note that the master integrals 



which produce singularities in the i5j~ sector can depend 
on Sei = 2p_ • pi raised to a non-integer power. Symbol- 
ically, 



2Re 



e 7 



= Bi+B2isei)-', (76) 



where B12 are functions that can be Taylor expanded 
around the Sei = limit. Since Sei ^ £,iVii the second 
term in the above equation provides additional 0{e) con- 
tributions to the exponents of singular variables. Obtain- 
ing those exponents correctly is crucial for constructing a 
valid expansion of the real-virtual corrections in inverse 
powers of e. Because there are two terms in Eq. ([75]) , we 
must introduce two functions i^i.2(^i, ^71) to parameterize 
the matrix element. We write 



(Oh 



Fi (6 ,m) + Ci"vrF2 (Ci ,Vi)^ F{Ci,r^i) 



(77) 



Away from the singular points ^1 = and 771 = 0, the 
function F{£^i, rji) is obtained by computing one-loop cor- 
rections to the radiative decay Z — >■ e"'"e~7 using stan- 
dard techniques. In the singular limits, we must distin- 
guish between the two contributions. Inserting Eq. (fTT]) 
into the phase space and performing the plus-distribution 
expansion for the ^1 and rji variables, we generate many 
terms with ^1 = and/or 771 = 0. We now discuss how 
to obtain Fi^2 at these singular points. 

First, we consider the case ^1 = 0, which corresponds 
to the photon 71 becoming soft. In that limit, both tree 
and one-loop QED amplitudes factorize into the products 
of the cikonal current and the corresponding amplitudes 
with the photon removed (s^ : 



A.(0.1)) 



p- • ei p+- ei 



M 



(0,1) 



(78) 



^P--Pi P+-P1. 

where ei is the photon polarization vector. Using this 
result, it is easy to find the limit of the function i^(^i, 771): 



lim 

?i=0 



^Im^i 



(0)* 

(0)* 



(79) 



Because no terms that behave like ^ appear in this 
limit, we conclude that 



f^i(0,7;i) = F(0,77i), ^^2(0,771)^0. 



(80) 



The next step is the calculation of the collinear 771 = 
limit. It is much more involved. The factorization in this 
case is given in terms of splitting amplitudes [lOj : 



M 
M 



(0) 

Z^e e+7i 

(1) 

2'— e+7i 



Split(?_ 



(1) 



(81) 



SphtW ,_^A^^ 



(0) 
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Hence, the one-loop amplitude factorizes into the one- 
loop splitting amplitude times the tree photon-less ampli- 
tude, and the tree splitting amplitude times the one-loop 
photon-less amplitude. The relevant splitting functions 
were computed in Rcf. [loj . They are given in terms of 
"standard matrix elements" . For an e* — s- ea7h splitting, 
the situation here, the QED splitting amplitudes are 



Sab 

Split*^^ = -2(r3(z)Split(°' -r4(z)Split 



(2) 



(82) 



where 



(83) 



and the two functions r3_4(z) parameterize loop contri- 
butions to the splitting functions. We must square the 
splitting functions and sum over the polarizations of the 
final-state particles. We find 



where the two functions ^3^4(2) are re-scaled versions of 
r3^4{z) in Ref. They are 



U f( \ Of ^ 



r4 = 7: 



1 e^l-Sne) 



2(l-e)(l-26) 



/2, 



where (5^ = in the FDH scheme [3^ and 6b = ^ in 
conventional dimensional regularization (CDR) . The two 
functions /i^2 can be found in Ref. [Toj . 

Finally, we comment on the construction of the expan- 
sion of the real-virtual contribution in plus distributions. 
The key point is that after the expansion is performed, we 
are able to get rid of -^1,2(^1, J?!) hi favor of -F(Ci,'7i)- for 
all values of Ci,f?i except the singular ones. At the sin- 
gular points, we have unambiguous expressions for i^i,2, 
as shown above. 



IV. REGULARIZATION SCHEMES 



t'"' X Split(°) 
Split X Split 



Split^"^ X Split^^^ ^ — Pe-y{e,z), 

Sab 

2 z{l + z) 



(84) 



Sab 1-2 



We can now use Eqs. ([81 



?7i — > limit of the matrix element: 



z(l + z) 
Pej[e, z)r3[z) -\ ri(z) 



1 - z 



where sie = 2p_ ■ pi. The two functions ''3,4(2:) are pro- 
portional to s^^ = (m|^77(l — z)) (see Ref. [13). Con- 



sequently, the first term in the right hand side of Eq. ((85)) 
contributes to the function -Fi(^i, 0), and the second term 
contributes to function ^2(^170). We find 



E-TUz 



where further simplifications are possible since z = 1— ^1. 
The limit of the function F2 is more complicated. It reads 



„ 2e 



-2Re 



Z— >e^e Z— >e^e ^ ^ 



if z{l + z) \ 

— -Pe7(e, z)r3{z) + — 7-4(2) , 

nz \ i-z J 



We note that most of the formulae presented in the 
previous sections do not make reference to a particular 
regularization scheme. They are valid independently of 
the scheme. Only the splitting functions, the tree-level 
Z — > e~^e~"f amplitude, and the one-loop Z 



, to derive the coUinear and Z 



e+e 7 



amplitudes change upon switching the 



scheme. As an illustration, consider the various con- 
tributions to the triple-collinear primary sector in Sec- 
tion III Al The double-soft contribution to the function 
Fi in Eq. (jSH) is given by the product of the square 
of eikonal currents and the tree-level matrix element for 
the Z — > e+e^ process. The eikonal current is scheme- 

(85) independent, while if we choose to work with physical 
four-dimensional polarizations of the Z-boson, the ma- 
trix element for Z e+e^ becomes scheme-independent 
as well. 

The real-virtual corrections proceed similarly. As ex- 
plained in Section IIIII we require the tree and one-loop 
matrix element for Z — > e~'"e~7, the one-loop matrix el- 
ement for Z — )• e+e^, and the tree- and one- loop split- 
ting functions for e — > e + 7. All of these objects are 
scheme-dependent, but this scheme-dependence is well- 
understood. In particular, the scheme-dependence of the 
splitting functions is given in Ref. [loj , while the scheme- 
dependence of the one-loop non-singular amplitudes can 
be found in Ref. [H. 

(86) Hence, it appears that within the framework discussed 
here any regularization scheme is allowed. The only non- 
trivial, scheme-dependent contribution at NNLO that 
needs to be computed explicitly is the two-loop vir- 
tual corrections. We emphasize that in this framework, 
no 0{e) terms of the double real emission amplitude 
Z e'''e~77 or real-virtual amplitude Z — ^ e~^e~j need 
to be known through higher orders in e, even in CDR. 
It seems at first glance that the 0{e) contribution to the 

(87) Z e~^e~"f amplitude is required, since it can hit a 1/e 
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pole when the other photon has become colhnear, lead- 
ing to a finite contribution. However, it has been sug- 
gested recently that this term cancels when the double- 
real and real-virtual corrections are summed (sBj'^. We 
have checked this statement by comparing the result from 
summing the double-real and real-virtual contributions 
in two different ways: with the full 0{e) contribution re- 
tained, and with the 0{e) term instead replaced by its 
coUinear limit. The sum of double- real and real- virtual 
is identical in these two cases, indicating that this term 
does indeed not contribution to the final result for the 
cross section. We note that simply dropping the 0{e) 
term would lead to a mismatch between the squared am- 
plitude and the approximation we use in the collinear 
limits, causing a divergence in the integration. It is non- 
trivial to track exactly how the 0{e) contribution cancels 
against similar terms in the collinear splitting functions, 
but since collinear limits are universal, the replacement 
that we advocate above appears to offer an easy, practical 
solution. 



V. NUMERICAL CHECKS 

To prove the utility of this method, we compute the 
contributions of double-real, real- virtual and virtual cor- 
rections to the decay rate of the Z-boson into leptons. 
For simplicity, we take the coupling of the Z-boson to 
leptons to be vector-like, and ignore all diagrams which 
contain photon vacuum-polarization contributions or its 
unitary cuts. We compare separately the double-real and 
real-virtual corrections to the four-and three-particle cuts 
of the vector- vector correlator, which we obtain using the 
optical theorem. We have presented these contributions 
separately because of the significant numerical cancella- 
tions between the double-real radiation, the real-virtual 
corrections, and the two-loop virtual terms that occur 
when summing them to obtain the total NNLO correc- 
tion to the decay rate. The three- and four- particle cuts 
of three-loop master integrals required for such compu- 
tation can be found in Ref. [s^. We present this com- 
parison in the CDR scheme. Note, however, that the 
polarization vectors of the .Z-boson are not continued to 
d-dimcnsions, making the tree decay rate Z — > e"'"e~ e- 
independcnt. We also set = 1- We write 



1 + 3^I + («A^^(2) 
4 TT Vtt/ 



where (J*-^' is further split into three contributions 

^(2) _ j(2) , ^(2) , ^(2) 

" ^ "rr ^ "rv ^ "vv 



(89) 



(90) 



^ The reason for this cancellation is a wcU-undcrstood indepen- 
dence of NLO QED corrections to Z e+e~7 on the regular- 
ization scheme. 



The two-loop virtual correction can trivially be ob- 
tained from the known result for the two-loop fermion 
form-factor. For this reason we do not present it here. 
From the analytic computation based on the optical the- 
orem we find 



(2) _ 0.5 1.5 

^ £4 + £3 



1.7246 14.074 



^(2) 

"rv 



1 



3 

73 



3.1794 22.. 



24.228; 



32.94. 



(91) 



These results can be compared directly to our computa- 
tions based on the soft and collinear limits of the relevant 
matrix elements. Before we present the corresponding 
results, we note one complication. There are interfer- 
ence contributions contained in Z e+e^e+e^ that can 
not be disregarded. Typically, these interference parts 
of the four-fermion final state correspond to certain cuts 
of non-planar diagrams and, hence, become part of our 
check. The four-fermion interference contribution only 
contains collinear singularities, and can be analyzed in 
the same way as the double-real emission contributions 
discussed in Section II. Because of the existence of this 
contribution, we split the double-real result into e+e~77 
and e~^e~e~^e~ final states. From our numerical calcula- 
tion, we obtain 



r(2),4e 

'rr. 



J? P 1 I o 



0.1799 



- 1.79, 



0.5 



^RR 



"rv - g4 



'f-hl^^'-^^ 22.6m, 



3.179 



22.84 



32.97(3). 



(92) 



The sum of the two double-real emission corrections 
agrees with Ec^. ([gi]) . as does the real- virtual contribu- 
tion, indicating the correctness of the numerical results. 



VI. CONCLUSIONS 

In this paper we have described in detail a subtrac- 
tion scheme which enables fully differential calculation at 
NNLO accuracy. By combining several ideas present in 
the literature, including an FKS partitioning of the final- 
state phase space and sector decomposition, the universal 
singular limits of amplitudes derived over a decade ago 
can finally be used to obtain actual physical cross sec- 
tions. Our ideas are explained using the simple test case 
of Z — e~^e~ as an example. We discussed how to par- 
tition the phase space based on the coUinear-singularity 
structure of the matrix element, and presented the ex- 
plicit phase-space parameterizations from which the soft 
and collinear singularities can be extracted as poles in e 
using sector decomposition. The treatment of the real- 
virtual corrections is described in a way that generalizes 
to more complicated processes. Numerical results that 
check our techniques were presented. We have chosen 
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to work in the CDR rcgularization scheme, although the 
presented framework remains vahd in other schemes such 
as FDH. It has been pointed out that difficuhies exist 
when extending FDH to NNLO [13 ■ Ahhough they can 
be fixed (ssj . witli our current understanding the use of 
CDR imposes no additional technical difficulties, as dis- 
cussed in Section ITVl 

One point we wish to emphasize about the result pre- 
sented here is its generalization to more complicated pro- 
cesses. As mentioned earlier, one of the problems with 
earlier sector-decomposition based approaches to NNLO 
calculations was the need to completely reconsider the 
phase space and extraction of singularities upon changing 
the process. In particular, if one knew the NNLO cor- 
rections to Z — > e+e~, but wanted to study the NNLO 
corrections to Z — >■ e+e^7. now would have to start from 
scratch. That is no longer the case for the framework 
described here. Differential Z decay serves as a building 
block for handling all final-state singularities, as we now 
describe. 

We will consider the real-radiation correction Z — 5- 
e+(p+) + e"(p_) +7(pi) -|-7(p2) + 7(P3), the most diffi- 
cult contribution, for definiteness. Introduce the follow- 
ing partition of phase space: 



is irrelevant; no soft or collinear singularities are asso- 
ciated with this momentum. We can simply reuse the 
NNLO results for Z e+e~ to obtain the corrections to 
Z — > e"'"e~7. Adding additional photons to the final state 
only increases the number of partitions required. In this 
sense, Z — > e~ serves as a building block for extract- 
ing final-state singularities from any process. While we 
have demonstrated this for only one partition, it follows 
similarly for the others. 

We are excited about the possible applications of these 
ideas to more phenomenologically interesting processes. 
We believe there is significant potential for applying these 
ideas to the calculation of 2 — >• 2 scattering processes at 
the LHC, and we look forward to their continued devel- 
opment. 
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(ij)G(l,2,3) 



Appendix: Splitting functions 



Here, ^i^jt allows and pj to be soft, but not any other 
particles; it also only allows collinear singularities when 
Pi,Pj,Pk are collinear. Sij_ki allows only pi and pj to be 
soft, and also allows only the collinear limits pi \\ pk and 
Pj \\ pi- D is the sum of all 6. It is simple to construct the 
appropriate d functions, as discussed in Ref. [l|. Consider 
the partition with (5i2,+ for concreteness. The contribu- 
tion of this real-radiation correction to the differential 
cross section is schematically 



We collect here the splitting functions that we employed 
in this computation, in the CDR rcgularization scheme. 
For the e — )• e -1- 7 splitting, we have 



Pe7(e, 2) = - (1 + ^) - <1 - 

For the e — ?• 67172 splitting, we find Q 



(96) 



da 

da, 



= / dLips, 



e^e 717273 



\M\^d{0-Oo) 



S12 



D 



(94) 



where O is an observable being studied and Ai is the ma- 
trix element. In this partition, there is no soft or collinear 
singularity associated with p^. We should therefore be 
able to use the phase-space parameterization and singu- 
larity extraction described in this paper, which handles 
the double- unresolved limit of photons pi and P2- To 
make this manifest, rewrite the phase space of Eq. ([M)) 
as 



dLips, 



eTe- 717273 



= ds+^i2[dp3][dp+-i2]dUps^ 
X S'-'^\pz - P3 - P+-12), 



(95) 



where S-| 12 is the invariant mass of all final-state par- 
ticles except the hard photon. The parameterization of 



momenta in dLipSg 



is chosen to be the same as 



in Section III Al The form of in this parameterization 



Pe^fi^, Z, Zi, Z2) 



2sieS2, 



■Pl(e, Ze, Zi,Z2) 



P2[e,Ze,Zi,Z2j 

Sle 



6 -(1-6) 



£2e 
Sle 



where the functions Pi 2 read 



D / \ I ^ ^ z^ 
Fl[€, Ze, Zi, Z2) = Ze, \ 6 



Z\Z2 



Z\Z2 



+ (1^2), 
(97) 



P2(e,Ze,2l,^2) = — ^ i^+e^(l+Ze) 



Z\Z2 

e{zl + ZIZ2 + zj)'^—^. 

Z\Z2 



(98) 
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